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Abstract 

The theory of superposition rules for sokitions of a Lie system of first-order differ- 
ential equations is extended to deal with analogous systems of second-order and the 
theory is illustrated with the very rich example of Ermakov- like systems. 



1 Introduction 



The characterization of non-autonomous systems of first-order differential 
equations 

-^ = Y\t,x), z = l,...,n, (1.1) 
at 

admitting a superposition rule is due to Lie and such a problem has be- 
ing receiving very much attention during the last thirty years because of its 
very important applications in physics [2]- [13]. The theory has recently been 
revisited from a more geometric approach in [11] where the role of the super- 
position function is played by an appropriate connection. The main point is 
that this new approach allows us to consider partial superposition of solutions 
as well and, furthermore, it also allows superposition of solutions of a given 
system in order to obtain solutions of a new system. Our aim here is to show 
how such a superposition rule may be understood from a geometric viewpoint 
in a very simple but interesting case, the so called Ermakov-Pinney system 
[T5|IT6] as well as for other generalizations of it. Such a system is made of 
second-order differential equations but the theory developed by Lie can eas- 
ily be adapted to deal with such SODE systems. We find in this way room 
for implicit nonlinear superposition rules in the terminology of [TTfTB] and 
the so called Ermakov-Lewis invariants [19] appear in a natural way as func- 
tions defining the foliation associated to the superposition rule. Moreover all 
reduction techniques developed for Lie systems [6IT2] are also valid in these 
cases. 
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2 Systems of differential equations admitting a superposition rule 



The superposition rule for solutions of (11.11) is determined by a function $ : 

X = $(x(i),...,a:(„);fci,...,/c„) , (2.2) 
such that the general solution can be written, for sufficiently small t, as 

x{t) = <l>(x(i)(t),...,X(^)(t);fci,...,/c„) , (2.3) 

where {x(a){t) \ a = 1, ...,m} is a fundamental set of particular solutions 
of the system (11. ip and k = (fci, . . . , A;„) is a set of n arbitrary constants 
associated with each particular solution. As a consequence of the Implicit 
Function Theorem, the function $(x(i), . . . , ■ ) : — > can be, at 
least locally around generic points, inverted, so we can write 

/c = ^(x(o),...,X(™)) (2.4) 

for a certain function : M"'(™+i) — > M". Hereafter in order to handle a 
short notation we start writing a;(o) instead of x. The function \E', also called 
superposition function, provides us with a foliation which is invariant under 
permutations of the (m + 1) variables. The functions $ and are related by: 

k = ^($(a;(i), . . . h, . . . , fc„),a;(i), . . . ,a;(m)) • (2.5) 

The fundamental property of the superposition function is that as 

/c = ^(x(o)(t),a;(i)(t),...,a;(^)(t)), (2.6) 

the function \E'(a;(o), . . . , a;(m)) is constant on any (m + l)-tuple of solutions of 
the system (11. II) . This implies that the 'diagonal prolongations' Y(t, X(o), . . . , X(m)) 
of the t-dependent vector field Y{t, x) = l^*(t, x) d/dx\ given by 

m 

Y{t, X(o), . . . , X(m)) = Yait, X(a)) , t G M , 

a=0 

where 

n Q 

yait,x^a)) = Y.ni^^ia))jr:^ (2.7) 

i=l ^-^{a) 

are t-dependent vector fields on M"(™+i) which are tangent to the level sets 
of \E', i.e. the components \E'* are constants of motion. The level sets of ^ 
corresponding to regular values define a n-codimensional foliation T on an 
open dense subset U C ]R"(™+^) and the family {Y{t)^ t E R} of vector fields 
in M"(™+i) consists of vector fields tangent to the leaves of this foliation. 

Remark that, as pointed out in \1A\, for each (a;(i), . . . , X(m)) ^ R"™ there is one 
point (x(o), . . . , X(^m)) on the level set J-'k of this foliation JF corresponding 
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to A; = {ki, ...,kn)e M", namely, . . . /c),X(i), . . . G JF^ (cf. 

(12. 5p ): then, the projection onto the last m factors 

pr : (x(o),X(i), . . . G M"(™+i) ^ . . . G M"™ 

induces diffeomorphisms on the leaves J^k of J^. Such a foliation gives us the 
superposition principle without referring to the function if we fix the point 
a;(o)(0) (i.e. we choose a. k = (fci, . . . , fc„)) and m solutions X(i)(t), . . . ,X(m)it), 
then X(o)(t) is the unique point in M" such that (x(o)(t), X(i)(t), . . . , X(m)(t)) 
belongs to the same leaf of J-' as (x(o)(0), X(i)(0), . . . , X(m)(0)). This means 
that it is only JF really matters for the superposition rule. 

Lie's main result [T] can be expressed as follows: 

Theorem: The system U.l]) on a differentiahle manifold N admits a super- 
position rule if and only if the t-dependent vector field Y{t^x) can he locally 
written in the form 

Y{t,x) = j2b^{t)X^{x) 

where the vector fields Xa, a = 1, . . . ,r, close on a r-dimensional real Lie 
algebra, i. e. there exist real numbers Cap ^ such that 

[X^,Xp] = j2c^l,^X,, Va,/? = l,...,r. (2.8) 

7=1 



The number m of solutions involved in the superposition rule for the Lie system 
defined by (12. 8p with generic ba{t) is the minimal k such that the diagonal 
prolongations of to A^'^ are linearly independent at (generically) 

each point: the only real numbers solution of the linear system 

r 

^CaX„(x(„)) = 0, a = l,...,/c 

a=l 

at a generic point (a;(i), . . . , X{k)) is the trivial solution Cq = 0, a = 1, . . . , m, for 
k = m, and there are nontrivial solutions for k < m. Then, the superposition 
function \E' is made up of t-independent constants of motion for the prolonged 
t-dependent vector field Y. We shall call them first integrals. 

A possible generalization consists on considering foliations which are not of 
codimension n, or even choosing different Lie systems with the same associ- 
ated Lie algebra for defining the prolonged vector field. These facts will be 
illustrated with several examples. 
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3 SODE Lie systems 



A system of second-order differential equations 



r{t,x,x), 



can be studied tfirougfi the system of first-order differential equations 



We call SODE Lie systems those for which X is a Lie system, i.e. it can 
be written as a linear combination with ^-dependent coefficients of vector 
fields closing a finite- dimensional real Lie algebra. There arc many interesting 
examples of such SODE Lie systems and next section is devoted to introduce 
some particular examples. 

3.1 The 1-dim harmonic oscillator with time- dependent frequency 
The equation of motion is x = —uj'^{t)x, with associated system 




with associated ^-dependent vector field 




d 



d 



X — V 



< 



V — —uj'^{t)x 



(3.9) 



and t-dependent vector field 



X = v--u:\t)x 



d_ 
dv 



which is a linear combination X — X2 — u!'^{t)Xi, with 





such that 



[X^,X2]^2Xs, [Xi,Xs] 



Xi , [X2 , Xs] — X', 



2 , 
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where Xs is the vector field given by 



1/ a d\ 

Therefore X defines a Lie system with associated Lie algebra s[(2, M). Actually, 
the vector fields Xi are the fundamental vector fields corresponding to the 
usual linear action and the basis of 5[(2,M) given by the following traceless 
2x2 real matrices 



ai 











(-io\ 
















10 , 




i ij 



(3.10) 



This system has no first integrals, i.e. there are not ^-independent constants 
of motion. 



3.2 The 2-dim isotropic harmonic oscillator with time- dependent frequency 
The system of equations of motion is 



with associated system 



xi — —cu'^{t)xi 

X2 = -Uj'^{t)x2 

Xi = Vi 
Vi — —U!^{t)Xi 

X2 = V2 
V2 = -U!'^{t)x2 

and t-dependent vector vector field 

d 2/ \ "9 d 2/ \ "9 

^ = ^^1-^ {t)Xl ^ h V2^ UJ [t)X2 -K- , 

OXi OVi 0X2 OV2 

which is a linear combination, X — X2 — cu'^{t)Xi, with 



^ d d 
Xi = Xi- h X2Tr- 

OVi OV2 



^ d d 

X2 = Vi- h V2- 



dxi dx2 



(3.11) 



such that 



[Xi, X2] — 2 X3 , [Xi, Xs] — —Xi , [X2, Xs] — X2 , 
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where the vector field X3 is defined by 

,,1/9 d d 

^3 = 77 Xi- Vi- h X2- 



dxi dvi dx 



2 




Once again X defines a Lie system with associated Lie algebra sl(2, M). How- 
ever, in the preceding case there is no constant of motion and we shall see that 
in this last one X1V2 — X2V1 is a constant of the motion. 



3.3 Milne-Pinney equation 



We call Milne-Pinney equation the second-order non-linear differential equa- 
tion [TBPI] : 

k 

X = -uj\t)x + — , (3.12) 
x-^ 

where k is a constant. It describes the time-evolution of an isotonic oscillator 
[2T|22] . i.e. an oscillator with inverse quadratic potential [23]. This oscillator 
shares with the harmonic one the property of having a period independent of 
the energy [21], i.e. they are isochronous systems, and in the quantum case 
they have a equispaced spectrum [25] . 



The corresponding system of first-order differential equations is 



X = V 



V = —uj'^(t)x + 



x-^ 



and the associated t-dependent vector field 

X — ^ ( "^i ) ^ \ 9 
dx \ x^ J dv 

This is a Lie system because it can be written as 

X = L2-uj\t)L^, 



where 

are such that 



Li = X—, 
ov 



k d d 
dv dx ' 



[Li, L2] = 2L3, [L3, L2] = -L2, [La, Li] = Li 



with 



L -- (x— - — 

^ 2 \ dx dv 
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i.e. they span a 3- dimensional real Lie algebra q isomorphic to s[(2,M). Ac- 
tually, one can show that they are the fundamental vector fields associated 
with the basis fl3.10p relative to the following action of the group SL{2, M) on 
a point (xo,fo) G M^: 

If the matrix A in SL{2, M) is given by 

(a A 

A = 

then (x, f)) = {x,v) with 



X = sign X 



k + [{/3v + ax){6v + -fx) + k{6(3/x^)f 

\ {6v + -fxy + k{6/xy 



(3.13) 



/ r ^2 kS'' I ^ X 



2 / a;2 



with sign(x) = 



3.4 Ermakov system 



Consider the system |26]|2 



X = Vx 
Vx = —UJ^{t)x 

y = vy 



V. 



y 



yi 



with associated t-dependent vector field 



^ d d 
ox oy 



which is a linear combination with time-dependent coefficients, X = ~uj'^{t)Xi- 
X2, of the vector fields 



^ d d 



d did 
X2 = Vx— + Vy— + 



' dx ^ dy y'^ dvy 
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This system is made up by two Lie systems, which correspond, respectively, 
to the examples of the 1-dimensional harmonic oscillator (13. 9p and the Milne 
equation (13.121) . both closing on a sl(2, M) algebra with X3 given by 



1 / a d 

-^3=0 V: 



2 \ dx "^dvr 



d_ 

' dy 



d 



3.5 Generalized Ermakov system 



y = -^9{y/x) -uj^{t)y 

y.i 

that for the choice f{u) = and g{u) = 1 reduces to the Ermakov system. 

This system of second-order equations can be written as one of first-order 
equations by doubling the number of degrees of freedom by introducing the 
new variables and Vy-. 



X = Vr 



-uj^{t)x + ^f{y/x) 



x-" 



y = u 



y 



Vy = -LJ^{t)y+—g{y/x) . 



Such system determines the integral curves of the vector field 

^(^-uj\t)x + ^f{y/x)^ l-+(-u;\t)y+^g{y/x) 



^ d d 

ox ov. 



d_ 



which can be written as a linear combination 

X = N2-uj'^{t)Ni 

where A^i and N2 are the vector fields 

d d d \ d d 1 d 

Ni = X— + y—, N2 = V:,— + -^f{y/x)— + Vy— + —g{y/x)—, 
ovx ovy ox x-^ ovx oy ovy 



Note that these vector fields generate a 3-dimensional real Lie algebra with a 



8 



third generator 

d d _ d \ 
^ 2 \ dx dvx ^ dy ^^dvyj 

In fact, as 

they generate a Lie algebra isomorphic to s[(2,M). Therefore the system is a 
Lie system. 

The fact that the systems of preceding examples are Lie systems with the 
same associated Lie algebra means that they can be solved simultaneously in 
the group 5*1/(2, M) by the equation (see e.g. [6,12j) 

gg~^ = uj^{t) ai - a2 , 

where ai and a.2 are given in fl3.10p . 



4 The superposition functions for these examples 



Consider first the 1-dimensional harmonic oscillator (13.91) . In order to look 
for a superposition rule we should consider a system like in (13. lip and check 
whether the vector fields Xi and X2 are linearly independent in a generic 
point. There are many points in which one can choose non trivial coefficients 
Ai,A2 and A3 such that Ai Xi + X2X2 + A3X3 vanishes in such a point. On 
the contrary, if we introduce another copy and obtain the system (13.111) . only 
from the vanishing of such vector in a point such that Xi f 2 — X2 f 1 = and 
X3 f 1 — Xi f2 = we cannot say that Ai = A2 = A3 = 0, therefore m = 2 and 
consequently there is a superposition rule involving two particular solutions. 

Note that such system (13. lip admits a first integral because the function F 
given by F(xi, X2, fi, ^2) is such that X2F = iff there exists a function 
F(^, f 1, V2) with ^ = X1V2 — X2V1, such that F{xi, X2, f 1, ^2) = -^(^, Vi, V2), and 
then from the second condition, 

dF dF 

XlT. ^X2T^ = 0, 

OVi OV2 

we obtain the first integral, which corresponds to the angular momentum, 
F(xi, X2, f 1, f2) = X1V2 — X2V1, which can be seen as a partial superposition 
rule. Actually, if Xi(t) is a solution of the first equation, then we obtain for 
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each real number k the first-order differential equation for the variable X2 

, , (iX2 ■ I \ 

Xi{t) —— = k + Xi[t)X2 , 

at 

from where X2 can be found to be given by 

X2 (t) =k'xi{t) + kxi it) f . (4. 14) 



In order to look for the superposition rule which does not involve quadratures 
we should consider three copies of the same oscillator, and the extended vector 
fields Xi and X2 given by 

d d d d d d 

Xi=Xi- hX2^ X2 = Vi- \-V2T, 

ovi 0V2 ov 0x1 0x2 ox 



We can determine the first integrals F as solutions of XiF = X2F = 0. The 
condition X2F = says that there exists a function F : ^ such that 
F{xi,X2,x,vi,V2,v) = F{^i,^2,vi,V2,v) with ^i(xi,X2,x,t;i,t;2,t;) = xvi-xiv 
and ^2{xi,X2, x, f 1, f2, v) = XV2 — X2V, and the condition XiF = transforms 
into 

dF dF dF 

— ^^2^ — yx— = 0, 

ovi 0V2 ov 
i.e. ^1 and ^2 are first integrals (Of course, ^ = X1V2 — X2V1 is also a first 
integral. They produce a superposition rule, because from 

XV2 — X2V = ki 
XiV — ViX = k2 

we obtain the expected superposition rule for two solutions: 

X = CiXi + C2X2 , V = CiVi + C2V2 , Ci = —, k = X1V2 - X2V1 . 



As a second example, consider the Milne system given in fl3.12p The generators 
of this Lie system with algebra sl{2, M) span a distribution of dimension two 
and there is no first integral of the motion for such subsystem. By adding 
the other s[(2,M) linear Lie system appearing in the Ermakov system, the 
harmonic oscillator with time dependent angular frequency, as the distribution 
in the 4-dimensional space is of rank three, there is an integral of motion. The 
first integral can be obtained from LiF = L2F = 0. But LiF means that 
F{x, y, Vx, Vy) = F{x, y,^) with ^ = xvy — yv^, and then L2F = is written 

dF dF X dF 
OX OX y-^ ot, 
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and we obtain the associated system of characteristics 
xdy — ydx dC, d{x/y) yd^ 



^ X ^ X 

from where the following first integral is found [T9] : 



+ ^=0, 



' x\'^ f x^ ^ 



^lj{x,y,v^,Vy) = +^ ^\y) ^^^^y^y^^y 

which is the well-known Lewis-Ermakov invariant [261I27II28I. 



We can follow a similar path in the case of the generalized Ermakov system. 
There exists a first integral for the motion, F : ^ M, for any a;^(t), because 
this Lie system has an associated integrable distribution of rank three and the 
manifold is 4-dimensional. 

This first integral F satisfies NiF = for i = 1, . . . ,3, but as [A^i, N2] = 2N3 
it is enough to impose NiF = N2F = 0. Then, if NiF = 0, 

OF dF 

dv^ dVy 

and according to the method of characteristics we obtain: 

dx dy dvx dvy 
X y 

and therefore there exists a function F : M'^ ^ M such that F{x,y,Vx,Vy) = 
F{x, y,^ = xVy — yVx)- The condition N2F = reads now 

dF OF ( y , , X , , \dF 



We can therefore consider the associated system of the characteristics: 

dx dy d^ 

Vx Vy -^f{y/x) + ^g{y/x) 

But using that 

—y dx + xdy dx dy 

^ Vx Vy' 

we arrive to 

—y dx + X dy (i^ 
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and integrating we obtain the following first-integral, with u = y/x, 

1 „ (l 



1 fU 



1 



dc = c. 



This first integral allows us to determine, by means of quadratures, a solution 
of one subsystem in terms of a solution of the other equation. 



4-1 The Pinney equation revisited 



We mentioned before the possibility of obtaining solutions of a given system 
from particular solutions of another related system. We next study a particular 
example, studied by Pinney long time ago [16]. Consider the system of first- 
order differential equations: 



which corresponds to the vector field 



X = Vx 

y = vy 



Z = V,, 



Vx 

Vy 

Vz 



-uj^{t)x + — 
x-^ 



X = v —+v — + v — 

^ dx ^ dy ^ dz ' x^ dv. 



d 



d 



d_ 



The vector field X can be expressed as X = N2 — uj'^{t)Ni where the vector 
fields A''i and N2 are: 



d d d 

dvy dvx ovz 



did d d 

N2 = Vy— + — h Vx^ + V^ — , 

oy X'^ oVx ox oz 



These vector fields generate a 3-dimensional real Lie algebra with the vector 
field A^3 given by 



N- 



d 



d 



d 



^ dx dvx dy dvy ' ~ dz ^ dvz j 

In fact, they generate a Lie algebra isomorphic to sl(2, M) because 
[iVi, N2] = 2iV3, [iVg, iVi] = iVi, [iVa, N2] = -N2 . 



d 



d 

+ z- Vz 



d 
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The distribution generated by these fundamental vector fields has rank three. 
Thus, as the manifold of the Lie system is of dimension six we obtain three 
time-independent integrals of motion. 

• The Ermakov invariant Ii of the subsystem involving variables x and y. 

• The Ermakov invariant I2 of the subsystem involving variables x and z 

• The Wronskian W of the subsystem involving variables y and z has 

They define a foliation with 3-dimensional leaves. We can use this foliation for 
obtaining in terms of them such a superposition rule. 

The Ermakov invariants read as: 

A = ^ (^{yVx - XVyf + k j , ^2 = ^ i^XV^ - ZV:,f + k ^ , 

and W is: 

W = yv:,- zvy . 

In terms of these three integrals we can obtain an explicit expression of x in 
terms oi y,z and the integrals Ii,l2,W: 



x = ^ [l2y' + hz' ± V4/1/2 - kW^ yz 



1/2 



This can be interpreted, as pointed out by Pinney [16j, as saying that there is 
a superposition rule allowing us to express the general solution of the Milne- 
Pinney equation in terms of two independent solutions of the corresponding 
harmonic oscillator with the same time-dependent angular frequency. 



5 The reduction technique for SODE Lie systems 



We aim to illustrate by simple examples the usefulness of the reduction tech- 
nique for Lie systems in this particular case of SODE systems. The main point 
is that given a Lie system in a homogeneous space for the associated group 
G, the knowledge of a particular solution allows us to reduce the problem to 
a new Lie system in the corresponding stability subgroup [61IT2] . Actually, if 
the Lie system is defined in the Lie group G by the equation 

g{t) g-\t) = ^{t) , ^7(0) = e (5.15) 

and we know a particular solution Xiit) of the associated system in a ho- 
mogeneous space, we can choose a curve g{t) in G such that ^(0) = e and 
<l>(^(t),xo) = x{t), and therefore, as also $(5'(t),xo) = x{t), there is a curve 
h{t) in Gx„ such that g{t) = g{t)h{t). Then h{t) is the solution of the Lie 
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system in Gx^ 

hit) h-\t) = Ad {g-\m^{t)) + r\t) m , /^(O) = e . 

Once such a Lie system has been solved the solution in the group is g{t) = 
g(t)h(t) and therefore the solution starting from any point yo of the associated 
Lie system in a homogeneous space is given by ^{g(t)h(t),yQ). As it is shown 
in the reduction technique, the transformation through the curve g{t) reduces 
the problem from the group G to the stability group of the initial condition 
of the particular solution xq. 

Instead, we can be interested in the reduction of the Lie system to one in the 
stability subgroup of a previously fixed point zq and we should proceed as 
follows. Let go be a fixed element in G such that ^{gojXo) = zq. Then g{t) 
is a solution of (15.151) iff g{t) = g{t) g^^ is a solution of the same equation 
as in fl5.15p but with ^(0) = g^^. Moreover, such curve g{t) can also be used 
to define the solution of the system in the homogeneous space starting from 
Xq by means of ^{g{t), zq). Note that if Xi{t) is a particular solution in the 
given homogeneous space with initial condition Xi(0) = Xq, we can choose a 
curve gi(t) in G such that ^{gi(t), Zq) = Xi(t) and then from ^{gi(t), Zq) = 
^{g(t),ZQ) = Xi(t) we obtain that: 

zo = <!>{g^\t),<l>{gi{t),zo)) = <l>{g^\t),^g{t), ^o)) = Hgi\t)g{t), zo). 

Thus, h(t) = gi^(t)g(t) lies in the stability group of zq, Gz^, and, consequently, 
gi{t) and g{t) differ on a curve in G^^. Recall that g{t) satisfies the equation 
(I5.15P but with ^(0) = g^^ and therefore h(t) is such that: 



h{t)h-' = j^{g^\t)m){9l\mt))-' 

= gi\t)gi{t)+Ad{gi\t)){s.{t)) = a'(t) G TjG,,, (5.16) 

and satisfies the initial condition h{0) = g^^{0)gQ^. 

Then, the transformation given by gi^(t) changes the initial equation in G as- 
sociated to a(t) into a new equation in Gz^ associated with a'(t) independently 
of the initial condition of the particular solution xi{t). When such stability 
group GzQ has a solvable Lie algebra g^^^ we can solve the equation in Gz^ for 
any a!(t). 

Conversely, if h{t) is a curve in Gzg solution of the Lie system in G^o 

hit)h'\t) = a'it), (5.17) 

satisfying the initial condition h{0) = gi^{0)gQ^, then g(t) = gi{t) h{t) go is the 
curve solution of our initial Lie system. On the other side, the equations defin- 
ing Lie systems are right-invariant and therefore if h'{t) is a curve solution of 
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f l5.17p but with the usual initial condition ^'(0) = e, then the solution of (15.171) 
with initial condition h{0) = gi^{0)go^ is h(t) = h'{t)gi^{0)gQ^ . Therefore, 
the solution of the initial Lie system is given by g{t) = gi{t)h'(t)gi^{0). 

Let us apply this theoretical development to three particular cases of this re- 
duction procedure: a reduction of a harmonic oscillator with a time-dependent 
frequency using one particular solution, a reduction of a Milney-Pinney equa- 
tion through a particular solution, and a reduction of a Milney-Pinney equa- 
tion by means of a particular solution of a harmonic oscillator with the same 
time-dependent frequency. 

Consider first the simple example (13. 9p of the harmonic oscillator with a time- 
dependent frequency and assume that a particular solution xi(t) is known. 
This is a Lie system in which is a homogeneous space for the group 
5'L(2,R). The orbit of the point (1,0) under the linear action is all R^, the 
stability group of such point being the 1-dimensional Lie subgroup generated 
by vd/dx. We can choose the matrix 



as a curve mapping the point (1,0) onto the given solution. Then we should 
make the change of coordinates corresponding to the transformation in the 
homogeneous manifold induced by gi^{t): 



and then we obtain the new second-order differential equation which is a first- 
order system in z, 



and allows us to find the general solution by means of a quadrature, because if 
u = z, the general solution of Xi{t) ii + 2 xi(t) u = is m = k/{xi{t)Y, for any 
constant and then we obtain by a second quadrature the expression (I4.14p . 
The fact that the new equation is solvable is related to the fact that the new 
equation is related with the solvability of g^o- This allows us to obtain its 
solutions by quadratures. Also, note that this change x = xi{t) z corresponds 
to the traditional d'Alembert method of reduction of order (see e.g. |34j). 

Consider now the example (I3.12p of the Milne-Pinney equation. This is a 
homogeneous case if we consider that either x > or x < and we restrict 
ourselves to one of these cases. Now, given a particular solution xi(t) of the 
Milne-Pinney equation the gi{t) constructed with xi{t) as in (15.180 transforms 
the point Zq = (1,0) into a solution {xi{t),Xi{t)) of this differential equation. 



xi{t) 

^Xi{t) Xi^{t) 




(5.18) 




xi{t) z + 2 xi{t) z = 
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Thus, as in the last case, wc can nse again gi^{t) to transform the initial 
equation in the group given by a(t) into a new one given by a'(t) that actually 
is a Lie equation in q^^. This new one will be solvable because 0^:0 solvable 
in this case. 



In this way, we start with: 



(5.19) 



The new equation in the group determined by the curve a!{t) constructed with 

hit)h-'it) = Ad{g^'itm4t)) + g-' g,{t) = ■^.'{t) , (5.20) 
turns out to be in this case: 



h{t) h-\t)^a!{t) 



xl{t) 



(k ai - 0,2) . 



If we consider a new variable r defined by 



T — 



lo xliO ' 

the new equation in the group is given by: 

dh{T 



dr 



-h-Ur)^k 



ai - a2 , 



(5.21) 



(5.22) 



then the solution with h{0) = e is /i(r) = exp((A;ai — a2)r) and we arrive to the 
solution for the original Lie system g{t) = gi{t)h{T)gi^{0). Thus, the solution 
for the time evolution is: 



xi{t) 



I 



cos(v^r) -^"mi 



—Vksm{y/kT) cos(\/fcr) 



a;r^(0) ^ 
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V 



Xi(t) COS (Vkr^ 



xi{t) sin(\/feT) 



VksiniVkT] / /t" \ • / \ cosl^/kr) sin( \/fcT)ii(t) 



-xi{0) xi{0) 

^ xi\0) ^ 
y-xi{0) xi{0)j 



Now, if $ denotes the action for the Milne-Pinney equation and we define: 




x^\0) 
-xi{0) xi(0) 




(5.23) 



then the set of solutions of the Milne-Pinney equation is written in terms of 
A, B as: 



xit) = 



' (52 + ^ + A'^k + {-B^ -k^ + A^) cos{2VkT{t)) + 2ABVksm{2VkT{t))xlit) 

2k 
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As a final example, consider once again the Milne-Pinney equation but assume 
that a particular solution xi{t) of the time-dependent frequency harmonic os- 
cillator with the same frequency as the Milne-Pinney equation is known. Re- 
call that it has been shown that both equations are related with the same 
Lie system in the same group G. In this case, once again gi{t) given by (15.181) 
transforms (1,0) into a particular solution {xi{t),xi(t)) of the time-dependent 
harmonic oscillator as a first-order differential equation. Thus, the transforma- 
tion induced in G changes the initial differential equation in G corresponding 
to the harmonic and the Pinney equation, which is the same one, into a new 
one in the stability group G^g. Now, the new equation in the group G^^ is 
given by: 



h{t)h-\t)= Ad {g^\t)){a{t))+g,{t)g^\t) 



1 



xlit) 



&2 



(5.24) 



and making use of a reparametrization given by (15.211) we obtain that the new 
equation is: 



dh{T) 

dr 

and the solution with h{0) = e is: 



h-Hr) 



-a2, 



(5.25) 



h{r) 



1 T 
1 



(5.26) 



from where we obtain the solution of our Lie system: 



git) 



xi{t) 

Xi{t) Xi^{t) 
Xi{t) Xi{t)T 



1 r W x^\0) 
1 ) \^-Xi(0) xi(0) 

^1/xiiO) 



^xi{t) xi{t)T + xY\t) ] \-i;i(0) xi(0) 

^ xi(t)/xi(0) -xi(t)xi(0)r xi(0)a;i(t)r 

^xi(t)/xi(0) - Xi{Q)/xi{t) - rxi(t)xi(0) xi(0)(xi(t)r + x^^it)) 

If we introduce the parameters A and B by 

/ 




xr^(o) 

-xi(0) xi(0) 




(5.27) 



where $ denotes the action for the Milne-Pinney equation, we obtain: 
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x{t) = + 2A^BT{t) + {A^B^ + k)T{ty 
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